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MONTE CARLO ANALYSIS OF LOBULAR GAS-SURFACE SCATTERING 
IN TUBES APPLIED TO THERMAL TRANSPIRATION 
by Jerry D. Smith and Charles A. Raquet 
Lewis Research Center 

SUMMARY 

A model of free molecular gas flow in tubes has been developed which combines a 
temperature dependent lobular reflection distribution (Logan-Keck-Stickney) with diffuse 
reflection at the wall, using Monte Carlo calculation techniques. The model is used to 
explain previously observed deviations in the thermal transpiration ratio from the 
Knudsen limiting law. These deviations suggest that when two reservoirs at tempera- 
tures T^ and T ^ are joined by a tube, gas molecules can be more likely to traverse 
the tube in the hot-to-cold direction than in the cold-to-hot direction; that is, Q 21 > 
when Tg > Tj, where is the tube transmission probability from reservoir i to 
reservoir j. 

The model yields magnitudes of transmission probability ratios Q 21 /Q 12 * or 
in Pyrex tubing in the range 1.09 to 1. 14, consistent with some experimental thermal 
transpiration ratios with T^ = 77.4 K, Tg = 295.0 K. The model also yields some ob- 
served effects of wall -temperature distribution, tube surface roughness, tube dimen- 
sions, gas temperature , and gas molecular mass. 


INTRODUCTION 

It has generally been assumed that the pressures measured in two regions of a 
vacuum system in the free molecular range are related to their associated temperatures 
by the following relation, which will be referred to as the Knudsen limiting law (refs. 1 
and 2) 
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However, several investigators (refs. 1 and 3 to 11) have noted deviations in the thermal 
transpiration ratio R = P^/Pg from this equation. (A list of symbols and their defini- 
tions is contained in the appendix.) 

When the two regions are joined by a tube, this law (eq. (1)) results from the follow- 
ing argument: Assuming that gas entering the tube is Maxwellian, the requirement that 



Figure 1. - Origin of conventional thermal transpiration ratio. Qjj is 
the probability that a gas molecule in chamber i ( which enters the 
connecting tube through plane Aj, will exit the tube through plane 
Aj (i, j = 1,2). 


the net mass flowthrough the tube be zero (see fig. 


1) results in the following equation: 


R = 




“ a T R m 



where Q^. is the probability that a gas molecule is chamber i which enters the tube 
through plane A^ will exit the tube into chamber j through plane A... (See fig. 1.) 

The quantity a T = C> 2 lA ^12 is the ratio of tube transmission probabilities in the two 
axial directions. If Qjg = Q 2 i> ^ en equation (2) yields equation (1). 

It can be shown that the transmission probabilities and Qgj are equal if it is 

assumed, for example, that the gas-surface interaction at the tube wall is any one of 
simple temperature independent types such as the following (ref. 12): * 

(1) Completely diffuse reflection (according to the cosine law) 

(2) Completely specular reflection j 

I 

(3) Some fixed degree of specular reflection regardless of the gas molecule's his- .1 

tory; the molecule has a constant probability of being specularly reflected each j 

time it strikes the wall; otherwise, it is diffusively reflected. This is j 

Maxwell's interpretation of gas-surface interactions (ref. 13). f 
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It has generally been assumed that for all gases, surfaces, and temperatures en- 
countered in ordinary vacuum technology, the reflections are completely diffuse. Con- 
sequently, equation (1) has often been used for calculating the thermal transpiration ratio 
(TTR). However, some experiments (refs. 1 and 3 to 11) indicate that under certain 
conditions a T > 1.00, implying that Q 21 >Q 12 > when T 2 T l - We wil1 refer to the 
experimentally measured quantity (P^/Pg)/!^ = R/ ^m as a E’ and the quantity 
Qgj/Qjg, which will be calculated theoretically in this report, as a T - In particular, 
Hobson, Edmonds, and Verreault found, in the most extensive experimental studies of 
thermal transpiration to date, that a E had values as large as approximately 1.31, for 
noble gases by Pyrex tubes, with = 77.4 K and Tg = 295.0 K (refs. 8 to 11). In 
order to explain such behavior, one must assume a temperature -dependent reflection 
pattern at the smooth tube wall (different from the three temperature -independent re- 
flection distributions described above), which results in a molecule having a greater 
probability of traversing the tube from the hot end to the cold end than in the opposite 
direction (ref. 9). 

Edmonds and Hobson suggested that different degrees of specular reflection of hot 
and cold molecules could account for the directional difference in Q^. (ref. 9). In line 
with this hypothesis, Miller and Buice (ref. 12) presented a temperature -dependent 
model for free molecular flow in a tube which explains the existence of the observed TTR 
deviations. Their model combined cosine law reflection (with complete temperature 
accommodation) with a degree of specular reflection (with no accommodation) dependent 
on the gas molecule energy. However, the transmission probabilities they calculated 
result in values of ranging from approximately 1.33 to 2.18 for tubes of L/r 
(length/radius) varying from 2 to 50, respectively. In general, these theoretical values 
(a^) were larger than experimental values (a E ), especially at higher L/r. Also, the 
trend of a T increasing with increasing L/r, does not agree with data from Hobson and 
his co-workers (refs. 8 to 11), which indicate that a E eventually levels off and prob- 
ably decreases, not increases, with increasing L/r. It should be noted that the ob- 
served TTR deviations from the Knudsen law depend on parameters other than L/r, 
such as tube surface, the wall -temperature axial distribution, and the geometry of the 
reservoirs behind the tube entrance planes (planes and A 2 in fig. 1). These pa- 
rameters were not extensively investigated by Miller and Buice. 

In another approach, Hobson (refs. 11 and 14) has shown that a pipe with unequal 
conductances in two directions can validly describe an accommodation pump. Briefly, 
this pump multiplies the pumping rate of a single tube with differential transmission 
probabilities by using many such tube sections with their ends maintained at cold and hot 
temperatures. By analyzing observed pressure changes with time, he derived values of 
the two different conductances through the entire pump. The ratio of these directional 
conductances, agrees well with the ratio of equilibrium pressures at the ends of the 
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pump. But, Hobson does not consider in detail the microscopic origins of this differ- 
ence in conductances (or, transmission probabilities) by dealing with the physical nature 
of the gas -surface interactions at the tube wall. 

It is the purpose of this report to investigate possible origins of the observed in- 
equality of transmission probabilities in single tubes with ends at different temperatures, 
using a gas-surface interaction model. A new model of rarefied gas flow in a cylindrical 
tube has been constructed. Gas molecule trajectories within the tube and transmission 
probabilities are calculated by Monte Carlo techniques. The basic difference between 
the approach presented herein and previous Monte Carlo studies is the following: Monte 
Carlo studies in the past have assumed complete diffuse reflection at the tube wall, or 
varying degrees of specular reflection (e.g. , refs. 15 to 19). In the present model, 
each gas molecule approaching the surface is assigned a certain probability of undergoing 
a diffuse reflection. If not reflected diffusely, it is re-emitted in a direction chosen 
from a lobular distribution function, the shape and location of which (relative to the spec- 
ular direction) is dependent on the gas and local surface temperatures. It will be shown 
that this temperature dependence results in the present model yielding a higher trans- 
mission probability in the hot-to-cold direction than the cold-to-hot direction. 


PRACTICAL IMPORTANCE OF THERMAL TRANSPIRATION 
EFFECTS AND THE PRESENT MODEL 

The importance of understanding TTR deviations from the Knudsen limiting law and 
their microscopic origins is twofold at the present time: 

(1) The common use of the low pressure thermal transpiration effect is to deduce the 
pressure in one part of a vacuum system, perhaps inaccessible to a pressure gage, 
knowing the pressure in another part and the temperatures in the two parts. Such cor- 
rections are often made throughout the pressure range (free molecular to continuum) 
using analytic descriptions of the thermal transpiration ratio that contain the experi- 
mental low pressure limit R m as a constant (e.g. , Liang's formula, discussed in 

ref. 9). If this experimental low pressure limit is not the conventionally assumed low 
pressure limit of R m = (Tj/Tg)^^, as is suggested in some experiments (refs. 1 and 
3 to 11), then analytic descriptions using R m yield incorrect values of the TTR at low 
(free molecular) pressures, as well as in the transition range. 

(2) The pumping action of a tube with differential conductances caused by gas-surface 
interactions, termed accommodation pumping by Hobson, is a new pumping principle at 
low pressures (ref. 11). It could be of practical importance, particularly for pumping 
gases of low molecular weight. 

The model in this report, utilizing the temperature dependent gas-surface reflection 
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distribution, presents another approach to the study of rarefied gases in tubes. Although 
the present study primarily deals with TTR effects, it may be also useful in interpreting 
other experiments involving gas -surface interactions. 

Consider, for example, a tubulated pressure gage calibrated in a facility where one 
would not expect temperature -dependent gas-surface scattering to occur in the tube. Such 
a situation might occur with relatively cold gases and a cold, microscopically rough 
metallic tube surface. If this gage were used in a facility where lobular scattering in 
the tube could occur, even from conventional engineering surfaces (refs. 20 and 21), 
such as within or near the hot (>1000 or 2000 K) gas of a thruster plume, the calibration 
may no longer accurately correlate gage reading with gas pressure or density. The 
pressure reading may be deceptively high because more gas molecules are entering the 
gas ionization chamber than entered the chamber when the calibration was done, with 
strictly cosine reflection occurring. 

In a related study, Ballance has considered various combinations of specular and 
diffuse reflections on the same molecular trajectory through a tube leading to the ioni- 
zation region of a pressure gage (ref. 22). This was done in an effort to explain differ- 
ences between neutral particle densities as determined from tubulated pressure gages 
and as determined from drag measurements on Explorer satellites containing the gages. 


PRELIMINARY ANALYSIS AND THE MODEL 

We will next describe in some detail the form of the lobular reflection pattern incor- 
porated in the present tube model. The use of this equation is the principal new feature 
in the model. While this lobe equation has had considerable success in predicting ob- 
served reflection distributions of molecular beams scattered from flat surfaces , it is a 
fairly recent (1966) addition to gas-surface theory and is probably not familiar to all 
readers. After this discussion, the tube model will be described in detail, showing how 
typical molecule initial positions on the tube entrance planes and initial trajectories are 
generated by Monte Carlo techniques and how molecules are traced until they exit the 
tube, being reflected from the wall either diffusely or lobularly. The number of mole- 
cules successfully traversing the tube is divided by the number of initial trials, yielding 
calculated transmission probabilities. 


Lobular Reflection Equation 

The lobular reflection pattern used in our model is given by the equation predicted 
by the Logan, Keck, and Stickney (LKS) hard-cube model (refs. 23 and 24): 
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where 

P(# r ) probability of reflection into d# r at 9^ 

T g ,T s gas, surface temperatures, respectively 

p nig/m g , ratio of masses of gas and surface atoms 

6 r angles of incidence and reflection measured from the surface normal, 
respectively 


and where 


B. = (- + ^ | sin 6 ■ cot 6 - ^ ~ ^ cos 9. 

1 V 2 / 1 r 2 1 


B„ = (— + M ) sin 9. esc 2 d 

2 V 2 / 1 r 


This equation predicts the variation of reflected flux intensity as a function of angle 
of reflection after gas molecules have interacted with surface atoms simulated by inde- 
pendent smooth cubes oscillating with a Maxwellian energy distribution. An important 
assumption, experimentally confirmed for many gas-surface combinations (ref. 24) , is 
that the tangential component of momentum of the gas molecule is unchanged during the 
interaction. 

Concerning the use of the LKS equation in the present tube model, it is important 
to note the following: Strictly speaking, equation (3) is valid only when the gas molecules 
approaching the surface are characterized by a Maxwellian energy distribution, with 
some temperature T . In the present tube model, each gas molecule, including those 

o 

reflected from the tube wall many times, is characterized by a temperature only, which 
is intermediate to the coldest and hottest wall temperatures. We do not deal explicitly 
with energy or velocity distributions, particularly of reflected molecules. To do so 
would greatly complicate the problem and is not necessary to explain several features of 
TTR deviations from the Knudsen limiting law. Stickney, Logan, Yamamoto, and Keck 
have shown that, based on the hard -cube model, little difference should be expected be- 
tween the scattering of a thermal Maxwellian beam and an appropriate monoenergetic 
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beam (ref. 25). Smith, Saltsburg, and Palmer have verified this experimentally, for 
argon and xenon beams on silver surfaces (ref. 26). They state that the dispersion in 
scattered beams in typical thermal energy gas-surface studies results not so much from 
the velocity distribution in the incident beam or from surface roughness as from the dy- 
namical nature of the surface atoms. The important experimental parameters are sur- 
face and gas temperatures. 

It is important to note that we use the gas temperature to indicate average gas mole- 
cular kinetic energy. Hence, when we refer to hot gas or cold gas, we will mean gas 
molecules of high or low average kinetic energies, respectively. We do not deal ex- 
plicitly with the spectrum of energies technically implied by a gas temperature. The 
relative lack of dependence of both observed and calculated reflected patterns on the 
velocity distribution in an incoming molecular beam partially supports our simplification 
of dealing strictly with temperatures. In the present model, we ascribe to each gas 
molecule the property of temperature, regarding temperature as a measure of average 
kinetic energy. 

There are several reasons for the use of equation (3) in the present study of thermal 
transpiration effects: 

(1) It contains explicit, experimentally verified dependences on gas molecular mass, 
surface molecular mass, angle of incidence, and gas and surface temperatures (see figs. 
2 to 4). The experimental results are from Hinchen and Foley (ref. 27), and the calcu- 
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Foley (ref. 27). 
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model (ref. 28). 


Figure 2. - Comparison of dependence of experimental and theoretical scattering distributions on surface temperature (from ref. 28). Argon on 
platinum; gas to surface atomic mass number ratio. 0.2; incidence angle. 67.^; gas temperature, 295 K. 
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(a) Experimental results for argon on platinum from Hinchen and Foley 
(ref. 27). 



(b) Corresponding theoretical patterns based on the LKS hard -cube model 
(ref. 28). 


: igure 3. - Comparison of dependence of experimental and theoretical scatter. ng distributions on the angle of incidence (from ref. 28). Argon on platinum; 
gas to surface atomic mass number ratio, 0.2; gas temperature, 295 K; surface temperature, 1081 K. 



(a) Experimental results for helium, neon, and argon on platinum 
(ref. 27). 



He J6CP 


(b) Corresponding theoretical results from the LKS hard-cube model 
(ref. 28). 


Figure 4. - Comparison of dependence of experimental and theoretical scattering distributions on mass of incident gas atoms (from ref. 28). Platinum 
surface; gas to surface atomic mass number ratio, 0.02 for helium, 0. 10 for neon, and 0.20 for argon; incidence angle, 67.5°; gas temperature, 

295 K; surface temperature, 1073 K. 



lated patterns are from simple hard-cube theory (ref. 28). These dependencies result 
in theoretical explanations of several aspects of TTR deviations. It should be noted that 
the theoretical distributions in these figures result from an earlier version of the hard- 
cube theory (nonclosed form, ref. 28), not from the use of equation (3). However, the 
scattering patterns given by this equation agree closely with those given by the analysis 
in reference 28 (ref. 23). 

Of particular importance in the present study are the recent experimental observa- 
tions of O'Keefe and Palmer of nondiffuse scattering of a 300 K helium beam from 300 K 
polished glass (ref. 29). For example, for 9 . = 70°, they observed a rather broad re- 
flection lobe with a maximum in the range 9 r = 40° to 60°. For comparison, a lobe 
maximum at about 62° is predicted by the LKS equation (eq. (3)), using p = 0.21 (to ap- 
proximate the helium-on-glass mass ratio), 9^ = 70°, and = T g /Tg = 1.00. (Their 
observed patterns also showed substantial reflected intensity at & r < 0°; i.e. , some 
molecules were backscattered, which is consistent with the assumption in the present 
model that some molecules are reflected according to the cosine law.) For argon on 
smooth glass, O'Keefe and Palmer noticed that the same nondiffuse trends, but to a 
much lesser degree. In general, for argon they corroborated results of Hurlbut 
(ref. 30), who observed predominantly diffuse scattering. 

(2) The LKS equation is probably the simplest closed-form reflection pattern other 
than diffuse reflection, which has experimental confirmation. This relative simplicity 
is quite important in terms of computer time because Monte Carlo simulation of many 
gas trajectories in very long tubes can be quite time consuming. 

(3) In general, cube models, which may seem grossly oversimplified to those famil- 
iar with the nature of real surfaces, have had considerable success in the treatment of 
gas-surface interactions. The hard-cube model, valid only for single gas-surface atom 
collisions and for the gas-to-surface mass ratios p less than 1/3, yields good scatter- 
ing distribution agreement for lighter mass gas species. A modification of this model by 
Logan and Keck involving the addition of a potential well (the soft-cube model, ref. 31) 
gives better scattering distribution agreement with heavier gas molecules (p > 1/3) and 
accounts for multiple collisions. However, the soft -cube equations are probably too 
complicated to use in Monte Carlo calculations. It might also be noted that the soft-cube 
model has been used to successfully predict energy accommodation coefficients for rare 
gases on tungsten (ref. 32). Also, Goodman has recently described a quantum mechanical 
basis for the cube models (ref. 33) . 


The Model 

The Monte Carlo approach to rarefied gas dynamics problems involves the genera- 
tion of many typical molecular trajectories according to certain assumed probability 


9 


distributions, and then appropriate averaging of the properties of the individual molecules 
resulting- from those trajectories and the various gas-gas (if any) or gas-surface inter- 
actions sustained. With gas-surface scattering dependent on temperatures and incident 
angles, whether an incoming gas molecule eventually exits the tube out the far end de- 
pends in a complicated way on its trajectory history. The calculation of transmission 
probabilities in the present problem is more amenable to a Monte Carlo technique than 
to integral equation solutions of the Clausing type . 

The basic idea in the present model of allowing molecules to be reflected at the tube 
wall either diffusely or lobularly is similar to Maxwell's useful intuitive concept of 
combining cosine and specular reflection (ref. 13). Trilling (ref. 34) has suggested 
combining cosine and hard-cube lobular patterns to describe some complicated observed 
patterns, arguing that a fraction of incoming molecules strike "hard-cube positions" on 
the surface while the remainder undergo multiple collisions either with a single surface 
layer atom or with microscopic roughness asperities and are hence re-emitted diffu- 
sely. Observed reflected distributions of molecular beams scattered from a small 
area of surface (such as the patterns for helium on polished glass, ref. 29) sometimes 
resemble superpositions of diffuse and lobular patterns. 

The calculation proceeds in a manner similar to other Monte Carlo treatments of 
gases in tubes. The basic theorem is the following (e.g. , ref. 35, p. 314): If a random 
variable | has the normalized probability density function f (u) , Uq < u < u^, then the 
probability density function of a random variable tj, defined as 



is uniform in the interval (0,1). Hence, if we wish to obtain a number belonging to a set 
of random numbers { S^} , which has the known or assumed density function f(u), we 
merely solve the following equation for the upper limit S^: 

/ f(u)du = R. (5) 

u 0 

where is a random variable with a uniform probability density function on (0,1). 

(As in other Monte Carlo treatments, values for R^ are provided by a computer sub- 
program in our calculations. ) 

We assume that the initial position of the gas molecule is uniformly distributed on 
the entrance plane of the tube and that the probability of being emitted in a particular 
direction from a surface (or, plane) is directly proportional to the cosine of the angle 
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(b) Step function wall temperature. 


I To = 295. 0 K 


z = 3L/4 


T x = 77.4 K 


(c) Partial ramp wall temperature. 



(d) Full ramp wall temperature. 

Figure 5. - Present tube model showing lobular reflection at the wall and 
three assumed wall temperature distributions. 


between that direction and the surface normal. Using appropriate forms for f(u) in 
equation (5) (uniform and cosine law probability density functions), it can be shown that 
one proper choice of initial coordinates and direction cosines (see fig. 5) is the follow- 
ing: 


= r ^R^ cos(27tR2) 

(6a) 

y^ = r^R^ sin(27rR2) 

(6b) 

z.^ = 0 

(6c) 

cos(A^) =^^2 cos(2ttR^) 

(6d) 
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cos(B^) = sin(27rR^) 


(6e) 


COS^) 


=yrri; 


(6f) 


With this initial position and direction, the intersection (x 2 ,y 2 , z 2 > °f the trajectory 
with the cylinder of radius r is calculated. If z 2 > L, then the molecule has exited the 
tube with no wall hits, and is counted both toward the direct tube transmission proba- 
bility Qp, the probability that a molecule passes through the tube with no wall hits, and 
the total tube transmission probability Q„, the probability that a molecule eventually 
leaves the tube through the exit plane, regardless of the number of wall collisions. If 
0<z 2 <L, then the simulated molecule has struck the tube wall, and decisions are then 
made as to whether the particle will be trapped permanently, or re -emitted diffusely or 
with the LKS lobe distribution, and what the re-emission temperature will be. Forthe 
present study, we have used a trapping probability of 0.0, since appreciable adsorption 
or surface diffusion is unlikely for the gases, temperatures, and tube sizes of present 
concern (ref. 11). 

A constant probability of a diffuse reflection F^ was assumed. The probability 
of a lobular reflection is then F L = 1 - F^. In this respect, the present model differs 
from that of Miller and Buice (ref. 12) and from the suggestion of Edmonds and Hobson 
(ref. 9) concerning hotter molecules having a greater probability of being specularly 
reflected than cold molecules. 

Following the gas-wall hit, the choice of a diffuse or lobular reflection is made by 
first choosing another random number, say Rg. If Rg < F^, we choose the reflection 
to be diffuse, determining the new direction cosines in the local coordinate system on the 
wall (with (x 2 , y 2 , Zg) as origin; see fig. 5) according to equations (6d) to (6f), only with 
newly chosen random numbers. Then, the intersection of this new trajectory with the 
cylinder is calculated. If Rg > F^, then the re -emission will be lobular, and the cal- 
culation of direction cosines is more complicated. 

The LKS probability distribution P(0 r ), equation (3) is a two-dimensional reflection 
pattern, in a plane formed by the incident trajectory and the surface normal at the point 
of impact. (See fig. 5.) To choose a reflection direction 6 r , O<0 r <77/2, measured 
from the surface normal in the appropriate plane according to P(# r ), the approach using 
equation (5) is not used. Although P(# r ) can be integrated in closed form, the integra- 
tion does not lead to a simple solution for the upper limit in equation (5). Instead of a 
numerical solution for this limit, the following procedure was used (based on, e.g. , 

ref. 35, p. 319): P(0_) was normalized by its calculated maximum value P(0 ), 

7 r ilia x 

yielding W; 0 < Pj^(0 r ) = P(9 r )/P(0 fflax ) < 1* Next, a pair of random numbers 
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was chosen, R^, R i+1# Letting 0 r = (7r/2)R i , we compare P N (0 r ) with R i+1 - If 
R^ + l < we accept the direction 0 r . The other two direction cosines in the local 

wall coordinate system follow from the known 0^ and the chosen 0 r , and the next cy- 
linder intersection can be calculated. If R^ + i > Pj<f(^ r ) » we reject the pair R^, R^ + ^, 
choose a new pair of random numbers, and repeat the process. We continue the process 
until we find an fL + ^ < P^($ r ) , or until 10 such attempts have been made at which time 
we abandon further choosing and merely send the molecule off at 0 rv , ov , its most prob- 
able value. It can be shown that such a procedure will statistically result in a choice of 
# r in accordance with the probability distribution function P(0 r ). It can be looked upon 
as ’’throwing a dart” at a unit square and accepting the coordinates so chosen only if the 
point lies below the curve P^(0 f ) (ref. 35, p. 319). 

When the molecule is reflected from the tube wall, its temperature is changed to 
Tqut> according to the following formulas using Knudsen’s familiar concept of the tem- 
perature accommodation coefficient (ref. 36): 


OUT “ a c|?s^ 

- t in] 

+ T IN’ 

for cosine reflection 

(7) 

OUT = a Lfy z) 

- t in] 

+ T IN’ 

for lobular reflection 

(8) 


where o/q and a-^ are constant temperature accommodation coefficients for cosine and 
lobular reflections, respectively, and T (z) is an assumed wall -temperature distribu- 
tion, like those shown in figures 5(b) and (c). For the initial gas-wall collision, t in 
will be either or Tg, depending on whether we are tracing molecules from the cold- 
to-hot and/or hot-to-cold end, respectively. Note that t out for one gas-wall inter- 
action will become t in for the next interaction, if there is one. 

After each new intersection is calculated, a test is made to see whether the molecule 
has left the tube through the entrance plane (z = 0), or out the exit plant (z = L), or has 
struck the tube wall again. If the molecule has left the tube, another initial molecular 
trajectory is generated according to equations (6a) to (6f), and the procedure is repeated. 

The entire process is repeated until a specified number of trials K q (typically 
10 000), have been completed. All trials are first performed for molecules entering the 
tube from the cold reservoir T^. Then, the entire process is repeated for K q addi- 
tional trials, for molecular entry from the hot reservoir Tg, with appropriate changes 
in initial gas temperature and T (z). 

The total number of molecules initially at T^ which eventually leave the tube 
through the exit plane is divided by K q , yielding a calculated value of Q^.. The total 
number of molecules leaving the tube directly from the entrance plane through the exit 
plane is divided by K q , yielding the quantity Qj-,. This direct transmission probability 
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Qp, not of particular interest in the present TTR studies, is independent of axial direc- 
tion because molecules directly transmitted never encounter the wall and, hence, never 
experience temperature -dependent reflection. (Therefore, results were used only 

as a check on the calculations and will not be presented.) The theoretical value of the 

1/2 

TTR, normalized by the temperature factor (T^/Tg) , is given by 


Q 


&ip — 


21 


Q 


12 


(9) 


Also calculated are the following: (1) the average number of gas-wall collisions per 
incoming gas molecule initially at T^, K^^k, j), which will be seen to also indicate the 
existence of TTR deviations from the Knudsen limiting law, (2) the axial distribution of 
these collisions, expressed as wall flux, (3) the axial distribution of average gas temper- 
atures of molecules striking the wall, and (4) the average temperatures of molecules 
leaving the tube through the two end planes. 

It should be noted that the calculated average temperatures are only simple averages: 
total temperature divided by total numbers of particles. Temperature, like any gas 
transport property, should properly be calculated using velocity weighting procedures. 
However, gas velocity magnitude changes at the surface and resultant reflected velocity 
distributions were not dealt with specifically in the present model. We have simplified 
the energy exchange by using temperature accommodation coefficients, as discussed in 
the next section. 


Use of Temperature Accommodation Coefficients 

The use of equations (7) and (8) involves a major simplification and departure from 
the LKS theory: simple, constant temperature accommodation coefficients dictating the 
reflected molecules' temperatures (or, energies). The LKS theory predicts a certain 
velocity distribution at any given reflected angle 0 , and, from this distribution, a most 
probable velocity at any angle 0 r (ref. 23). However, the LKS hard-cube velocity pre- 
dictions have not as yet had the success experimentally that the directional predictions, 
such as those shown in figures 2 to 4, have had. Some experimental determinations of 
average reflected molecule energies (i.e. , temperatures) indicate a variation of average 
reflected gas temperature with 0 with some trends different from LKS predictions 
(refs. 21, 37, and 38). Other reflected gas temperatures have been measured to be 
fairly independent of 0 r (e.g. , refs. 20 and 21). To the authors' knowledge, no reflec- 
ted gas energy distributions have been measured for light mass gases (such as helium) 
on glass surfaces. Only a few measurements of accommodation coefficients for light 
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mass gases on glass have been made (ref. 39). For most of the present work, we used 
the value a = 0.35, from reference 39. 

It should be mentioned that at present it is not unusual for a gas-surface interaction 
model to predict only certain gas-surface properties successfully. It is generally known 
that the cube models predict scattering patterns well and that some lattice theories (in 
which incoming gas atoms interact with a one-, two-, or three-dimensional lattice of 
atoms in the solid) successfully predict thermal accommodation coefficients (ref. 40). 
The interaction assumed in the present tube model can be regarded as a reasonable and 
convenient hybrid. We incorporate the LKS scattering direction predictions with the 
more conventional simpler ideas of cosine reflection and temperature change through 
accommodation coefficients, to yield a tractible Monte Carlo problem. Although the 
present model is convenient, it should be noted that it is not really self-consistent: the 
use of a temperature accommodation coefficient is inconsistent with the assumption in 
the LKS model that the tangential component of momentum is unchanged during the inter- 
action. However, even if the LKS velocity distributions were verified, we feel that 
sampling from both directional and velocity magnitude distributions would require pro- 
hibitive computer time for the present study. Fortunately, our calculated values of 
are relatively insensitive to a q and a-^, which could indicate that values of long tube 
transmission probabilities (hence, of a^,) are not very sensitive to gas temperatore (or, 
energy) changes at each gas -surface interaction. 


Statistical Accuracy of Transmission Probabilities 


As a measure of the statistical accuracy of Q„, 
or limit for is calculated 


the 95-percent confidence interval 


e N 


1.96 




(10a) 


This interval is interpreted as follows: Monte Carlo calculations of based on K q 
trials should be within of the exact theoretical answer 95 percent of the times the 
calculation is done, assuming that the answers are normally distributed about the 
theoretical answer. This is usually a reasonable assumption if is not extremely 
small. The theoretical value of should be used in equation (10a), but since its 
value for the present problem is unknown, the calculated value will be used. The most 
general estimate of Monte Carlo error (when no assumptions are made about the distri- 
bution of calculated Q^'s) i s given by the Chebyshev inequality 
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| (Monte Carlo) -Q^ (theoretical) | < 


(theoretical) Jj. - (^.(theoretical^ 

pk: 


= dr 


(10b) 


where p is the probability of falsity of the inequality in equation (10b). (From ref. 35, 
p. 12.) The quantity (1 - p) x 100 percent is referred to as confidence, and the quanti- 
ties and e r are confidence limits or intervals. For 95-percent confidence 

(p - 0.3b, 


*c = 4 - 47 


/q ! .( theoretical) D - Q^.(theoreticaI)| 


K 


(10c) 


Our use of e^(<e£,) I s 1° some extent justified by comparison of our results for trans- 
mission probabilities with previous analytical results, for 100 -percent diffuse reflection 
and for partial specular reflection at the tube wall. This will be discussed in the first 
part of the RESULTS AND DISCUSSION section, which follows. The important feature 
of all estimates of Monte Carlo error is the inverse dependence on the square root of 
the number of trials. Quadrupling the number of trials only halves the confidence in- 
terval. 


RESULTS AND DISCUSSION 

In the following sections, we first compare results of the present calculation of 
transmission probabilities with some previous theoretical results, where applicable, to 
show the accuracy of the present model. Then a brief summary of some previous ther- 
mal transpiration ratio (TTR) experiments is given, listing the specific parametric ef- 
fects which a theoretical tube model must attempt to explain. Next, preliminary to the 
results of the generation of molecular trajectories in the tube model, some properties 
of LKS lobular reflections are discussed, with emphasis on gas and surface temperature 
dependence. Then, we present results for tube transmission probabilities and results 
concerning numbers of gas-wall collisions, when both cosine law and lobular reflections 
occur at the tube wall. Finally, we show to what extent the model is capable of explain- 
ing the various experimental TTR effects observed. 
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Comparison with Previous Theoretical Tube Transmission Probability Results 


The results of the present model can be compared with previous studies for which 
all reflections from the wall are diffuse by setting = 1.0. Our tube transmission 
probabilities in table I, presented as functions of L/r, are within of all but one of 
other analytical and Monte Carlo results shown (refs. 41 to 45), particularly those of 


TABLE I. - TUBE TRANSMISSION PROBABILITIES FOR COMPLETELY DIFFUSE (COSINE LAW) 


REFLECTIONS AT THE WALL (F c - 1.0, IN THE PRESENT MODEL) 


Tube length 

Number of 

Tube trans- 

95-percent 


Previous values of Q.. 


to radius 

Monte Carlo 

mission 

confidence 






ratio. 

trials, 

probability 3 

interval, 

Richley 

DeMarcus 

Clausing 

Ward and 

Miller 

L/r 

K 

(present model), 


and 

(ref. 42) 

(ref. 43) 

Fraser 

(ref. 45) 


0 

Q.. 


Reynolds 



(Monte Carlo, 




13 


(ref. 41) 



with K q 









400 000) 









(ref. 44) 


0. 5 

40 000 

0.801275 

0.0039 

0. 801 

0.80127 

0.8013 


0. 80128 

1.0 

10 000 

.6772 

.0092 

.671 

.67198 

.6720 

0.6724 

.67209 

2.0 | 

10 000 

.5102 

.0098 

.513 

.51423 

.5136 


. 51472 

4.0 

10 000 

.3580 

.0094 

.355 

.35658 

.3589 

.3570 

.35629 

10.0 

20 000 

.19025 

.0054 


.19099 

.1973 

.1914 

.18656 


a Qj 2 Q 2 I r Qjj independent of direction when F^ - 1.0. 


TABLE II. - TUBE TRANSMISSION PROBABILITIES FOR 
PARTIAL SPECULAR REFLECTION AT THE WALL 


[Tube length to radius ratio, 2.0; number of Monte 
Carlo trials, 20 000; probability of specular 
reflection, 1.0 - Fj-,.] 


Probability of 

Tube transmission probability , a 

9 5 -per cent 

cosine law 
reflection, 

F c 


Q. . 
^13 

confidence 

interval, 

e N 

DeMarcus 
(ref. 47) 

Present model 

1.0 

0.51423 

0. 51625 

0.0069 

.8 

.58247 

. 58235 

.0068 

.6 

.65890 

.6588 

.0066 

.4 

.74693 

.7476 

.0060 

.2 

.85412 

.8535 

.0049 

0 

1.00000 

1.0000 

0 


a Q 12 = Qgj - Q^j, independent of direction. 
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DeMarcus, which are regarded as probably the most accurate (refs. 44 and 46). (The 
exception is Clausing's value for L/r = 10.0. But Clausing's values for higher L/r are 
believed now to be relatively inaccurate (ref. 44).) 

We can also essentially duplicate DeMarcus* results for partial specular reflection 
(ref. 47) by replacing the procedure of choosing from P($ r ) when a lobular reflec- 
tion is chosen, with the simple specular condition 9 r = 0.. Our results for L/r = 2.0 
are shown in table II as a function of F^, and are again within of DeMarcus' 
answers. 

Another comparison with previous theory concerns the use of an assumed constant 
nonzero trapping probability (or, sticking coefficient) at each gas-tube wall collision. 
Smith and Lewin (ref. 17) calculated tube transmission probabilities as a function of 
L/r and the sticking coefficient at the wall, using cosine law reflection for particles not 
trapped at the wall. The present model essentially reproduces these results also. (How- 
ever, as mentioned previously, for the present TTR studies, a sticking coefficient of 
zero was used.) 


Summary of Some Previous Experiments 


There are many factors that could affect TTR deviations. Several of these factors 





To pressure gage and 
gas addition system 



J Liquid level 
y height 


Liquid nitrogen 


\a i Two-bulb, or buib-bulb apparatus. 


lb) Bulb-elbow apparatus. 


— ► To gas addition system 




Leached Pyrex 
section — 


Dewa r 
motion 


‘V Smooth Pyrex section 


Liq u id 
nitrogen 


Liq u -d level 
height 


(c) The accommodation pump; three stages shown (elbow-elbow). 


Figure 6. - Schematics of types of apparatus used to investigate the thermal transpiration 
ratio by Hobson's group (refs. 8 to ID 
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have been determined experimentally, particularly in the extensive studies of Hobson's 
group (refs. 8 to 11 and 14). In the following brief summary of the effects suggested by 
these experiments, the TTR effects will be discussed in terms of the quantity 
a„ = R/R , which will be referred to as the temperature normalized thermal transpira- 

hj m 

tion ratio (experimental) . 

Effect of reservoir geometry . - Figure 6 shows the three types of apparatus used by 
Hobson's group to investigate the thermal transpiration ratio. These types are charac- 
terized by the "volumes" joined by the tubes being either glass bulbs, or 90 or 180° 
elbows. In general, the lower part of each apparatus was dipped into liquid nitrogen to 
varying levels, probably giving rise to a fairly sharp wall -temperature break in the 
tubes. The pressure in the cold region (P^) was derived from the pressure in the hot 
region (Pg) and some other measurable parameters assuming conservation of gas phase 
molecules (ref. 9). Generally, values of a^, in the elbow apparatus (figs. 6(b) and (c)) 
are larger than and more sensitive to liquid nitrogen level variations than a^ values in 
the two-bulb apparatus (fig. 6(a)) (ref. 9). For example, for helium with the liquid ni- 
trogen level half way up the tube, elbow apparatus values range approximately from 1. 13 
to 1.31, while two-bulb apparatus values are roughtly 1.03 to 1.11. 

Effect of wall -temperature distribution . - As stated before, for elbow apparatus, 
ag was sensitive to liquid -nitrogen level, increasing with increasing nitrogen height 
and then either leveling off as the height increased or going through a maximum (fig. 7 
and ref. 9). Note that, in general, a^, is still not equal to 1.00 even when the entire 
tube is at a single temperature, particularly when the tube is at 295.0 K (liquid -nitrogen 
height at 0 cm in fig. 7). For the two-bulb apparatus data, a^ is more insensitive to the 
liquid-nitrogen level (fig. 7 and ref. 9). The addition of a heater wrapped around the tube 
right above the liquid -nitrogen level, which probably sharpened the wall -temperature 
break, either slightly lowered the value of R (or, a^,) or left it unchanged (fig. 7, sam- 
ple 1 , with and without heater) . 

Effect of gas type . - All other factors held constant, for the noble gases, a^ de- 
creases with increasing gas mass. For a tube with L/r ^ 13.6, a^ varied from about 
1. 18±0.01 for helium (m^ = 4) to about 1.05±0.01 for xenon (m^ = 130) (ref. 11). 

Effect of tube material . - The results of Hobson (ref. 11) as well as those of Lund 
and Berman (ref. 48), suggest that metal tubulation, no matter how smooth or rough the 
surface may be, results in values for a-g within a few percent of the expected value 
1.00; tube transmission probabilities do not appear to be directionally dependent. So, 
the larger TTR deviations from Knudsen's limiting law seem at present to be restricted 
to smooth Pyrex tubulation. 

Effect of surface preparation . - Hobson found that leaching the Pyrex tubes (i.e. , 
microscopically roughening their surfaces by the action of HC1) caused the TTR devia- 
tions to virtually vanish. The greater the degree of leaching, the closer the TTR ap- 
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Thermal transpiration ratio (temperature normalized); 
a^ (experimental) and a T (theoretical) 



(a) Cosine law, or diffuse. 



Surface to gas 



(c) Theoretical LKS lobular distributions (eq. (3)). Gas to surface atomic mass 
ratio, 0.21; incidence angle, 60°. (T R = T$/Tg. ) 


Figure 8. - Three types of gas-surface reflection patterns. 
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r, 

L, 
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cm 

cm 
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45 
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45 
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Sample 1, without heater, 

Bulb -el bow 
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45 

47.4 

Sample 1, with heater 

Bulb-bulb 

.95 

24 

25.3 

Sample 1, without heater, 





after vacuum bake at 





500 C for 12 hr 

Bulb-bulb 

.95 

24 

25.3 

Sample 2 

O Theory 

1.71 

24 

14.0 

Some calculated values of 


(assumed) 

(assumed) 


a T s Q2l/3 12. from p re- 


sent model; see table IV 
for model parameters 


1.0Q I I 1 1 1 

0 5 10 15 20 25 

Liquid -nitrogen height in tube, cm 


Figure 7. - Normalized thermal transpiration ratio for helium as function of liquid -nitrogen height in tube. Data for curves 
from Edmonds and Hobson (ref. 9); Tj = 77.4 K; ^2 = 295.0 K; P^- 4xl0'6 torr. 


proached the anticipated theoretical value R m (ref. 10). 

Effe ct of preadsorbed gas layers . - When layers of Ar, Kr, or Xe were preadsorbed 
on the Pyrex tube, the value of a E for helium was unchanged (ref. 11). 

Effect of tube L/r . - This effect, although not clearly established, seems to indi- 
cate that the Knudsen limiting law (eq. (1)) holds for apertures (L/r = 0.0) and that, be- 
yond a certain L/r, a E decreases -with increasing L/r, perhaps back to 1. 00. The 
existence of a maximum (or of maxima) is hence suggested at some intermediate L/r. 

Effect of number of stages in a multistage pump . - Using the finding that leaching 
removes the directional difference in tube conductance, Hobson combined leached tubes 
with smooth tubes in the arrangement shown in figure 6(c), each leached-smooth tube 
combination forming a stage. The effect of increasing the number of stages was to in- 
crease the pumping action of a single smooth tube with the ends at different tempera- 
tures. For helium the value of the pump ratio P A /P B in the three-stage pump was 
2.10; in the 28-stage pump, 23.3. For a single-stage pump, the ratio, which is then 
theoretically equivalent to the previously defined a E (ref. 11), was 1. 18±0.01. 

Effect of hot chamber temperature Tg. - The warm half of a three -stage pump was 

changed from 295 to 600 K, resulting in an increase in pump ratio for helium from 2. 10 
to 2. 50. Although this increase cannot be reliably extrapolated to a single -stage pump 
or just one smooth tube (ref. 11), the qualitative trend of a E increasing with the hot- 
gas - hot -tube temperature was established. 


Lobular Reflection Distributions 

Shown in figure 8 are various types of possible reflection patterns: cosine law, 
specular, and lobular. In such diagrams as figures 8(a) and (c), the probability of emis- 
sion in a certain direction is proportional to the length of the radius vector in that direc- 
tion, as shown in figure 8(a). The LKS distributions (fig. 8(c)) depend on temperature 
only through the ratio T„ = T /T and on the nature of the gas and surface only through 

xv S g 

the mass ratio p = m /m . In figure 8(c) and throughout the paper, the value of 

g s 

fi =0.21 was used for helium on Pyrex (sodium borosilicate) , after estimating the av- 
erage mass of an atom in the glass surface at m =19.2 by weighting the masses of 
species within the Pyrex in proportion to their occurrence in the bulk. Considering 
the complicated structure of glass, p must be regarded to a certain extent as an adjust- 
able parameter. However, our results for a,p = Q21 //< ^12 are no * extremely sensitive 
to this parameter. 

It should be noted that the exceptionally good agreement between theoretical and ex- 
perimental reflection distributions in figures 2 to 4 occurs, coincidentally, in a situation 
where the mass ratio p = 0. 2 (argon on platinum). Since the LKS model depends on the 
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ratio of masses, it does not matter theoretically whether the ratio n = 0.2 arises from 

argon on platinum, or from helium on Pyrex with an assumed m = 19. 2. 

s 

Because the experimental arrangements being considered should result in a fairly 
sharp wall -temperature break, one would expect many collisions in the tube with 
T R = T g /Tg = 77. 4/295. 0 (hot gas on a cold wall) and with T R = 295. 0/77. 4 (cold gas on 
a hot wall). Figure 8(c) shows that the lobe for the hot gas on a cold wall is directed 
away from the specular direction toward the surface tangent and the lobe for the cold gas 
on a hot wall is directed more toward the surface normal. It is this differential dis- 
placement of lobular patterns, on either side of the specular direction, that results in the 
difference in tube transmission probabilities suggested by a R > 1.0. One would also 
expect gas -wall collisions in the tube with T R = 1.0 = 77.4/77.4 = 295.0/295.0, where 
the gas and wall temperatures are the same. However, there is no differential in lobe 
position for a hot gas on a hot wall and a cold gas on a cold wall, because the tempera- 
ture ratio is 1.0 in either case. 

Figure 9 shows the difference between the calculated angular position of lobular 
maximum 6 and the specular direction = 0. as a function of 0. for the 

HlaX. bp6C 1 1 

three values of T R mentioned previously. It is seen that the differential displacement 
is always present for 6 ^ > 0° and that the effect increases as the angle of incidence 
increases. In particular, for the cold gas on a hot wall (T R = 295. 0/77. 4) the 
lobes are displaced strongly toward the surface normal for large angles of incidence 
(large negative values of A 0 in fig. 9). For T R = 1.0, the lobe maximums are also 
displaced above the specular direction, the effect increasing with increasing 0^, but to 
a much lesser extent than lobes with T R = 295.0/77.4. 



Figure 9. - Difference between lobe maximum and specular direction (A0 = 0 max - 0 S p ec 
where 0 S p ec = 0j) as function of incidence angle (from eq. (3)). Gas to surface 
atomic mass number ratio, 0.21. 
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These properties of lobular reflections will be used in interpreting results for 
partial lobular reflection in tubes. 


Partial Lobular Reflection in Tubes - Transmission 
Probabilities and Gas-Wall Collisions 

The parameter , the probability of a gas molecule collision with the wall being 
followed by a diffuse reflection, must be regarded as an adjustable parameter, con- 
sidering what is now known about gas -surface interactions. Figure 10 shows the effect 
of varying F^ on the two transmission probabilities and For these calcu- 

lations, L/r = 14.0, which is typical of some tubes used by Hobson, Edmonds, and 
Verreault (refs. 8 and 11); p = 0.21, as discussed previously, K q = 10 000, which, by 
equation (10a), leads to the 95-percent confidence limits shown; chq = = 0.35, an ex- 

perimental value for helium on glass measured between about 250 and 400 K and rela- 
tively temperature independent in that range (ref. 39); and T (z) is the step function 
shown in figure 5(b). Also, we have assumed that F^, is the same constant for all gas- 
wall collisions, independent of gas or surface temperatures or angles of incidence. Note 


o e <= 



Figure 10. - Dependences of transmission probabilities, the ratio of transmission probabilities, and the average 
number of gas-wall collisions per incoming molecule on the probability of cosine reflection at wall. Tube 
length to radius ratio, 14.0; number of Monte Carlo trials, 10 000; temperature accomodation coefficient, 
o = o c = a L = 0.35; gas to surface atomic mass number ratio, 0.21; step function wail temperature distri- 
bution (see fig. 5(b)). 


23 



in figure 10 that, when F^, = 0.0 and all reflections are lobular, Q 12 = Q 21 = 1-000, as 
expected, since with no backscattering possible all molecules must traverse the tube 
starting from either end. Note also that when F c = l.o and all reflections are diffuse, 
Q 12- Q 21 a0 ' 147 ' since no temperature -dependent scattering has then occurred. Be- 
tween these two values of F c we find Q 21 >Q 12 and a T > 1-00, which is consistent 
with the experimentally observed TTR effect, that is, a^ > 1.00. 

The maximum value of a T in figure 10 is a T = 1. 119 when F c ^ 0. 3. In the 
range 0. 1 < F^ < 0. 6, values of a T are in the approximate range 1. 05 < a T < 1. 119, 
which is consistent with values from some two-bulb data (fig. 7), but somewhat less 
than the magnitudes of a^. from elbow data. We might expect our theoretical results 
to compare better with two-bulb data because molecular entry from a large bulb into the 
tube would be expected more Maxwellian, as assumed in our model, then entry through 
an elbow into the tube. 

As shown in figure 10, ^^(1,2) the average number of gas-wall collisions per in- 
coming molecule in the Tj - T 2 direction, is greater than Kq W ( 2,1). Thus, initially 
hot molecules T 2 make slightly fewer wall collisions, on the average, than the initially 
cold molecules Tj starting at the opposite tube end. Through the remainder of this 
report, the subscript 2 will refer to the hot temperature (gas or Avail), and the sub- 
script 1 Avill refer to the cold temperature. We have extrapolated Avith dashed lines in 
figure 10 the curves for K Q W (i, j) to their calculated common value of about 14.0 at 
F^ = 1.0 and to their approximate expected values at F^-, = 0.0. Note that when 
F^ = 0.0, even though all molecules ultimately transverse the tube starting from either 
reservoir, we would not expect the average number of gas-Avall collisions in the cold-to- 
hot direction K^^(l , 2) to equal K^^(2,l). The reason for this is the differential dis- 
placement of typical lobular patterns in the two different axial directions caused by dif- 
ferent typical wall -to -gas temperature ratios in those directions. The K^ w (i, j) curves 
go through a maximum and decrease to a common value of 14.0, as F^, goes to 1.0. 

As the probability of diffuse reflection increases, more molecules are backscattered out 
the entrance plane soon after they enter the tube and hence make fewer total gas -wall 
collisions. 

Figure 11 shows how the Avail flux varies axially, when = 0.3 and L/r = 14. 0. The 
wall -temperature break between 77. 4 and 295 K is located at L/2. The tube Avas divided 
into 35 equal axial segments. For 10 000 initial trials from the entrance plane, simu- 
lating an inlet flux of molecules entering from the 77. 4 K reservoir, the total number 
gas-wall hits in each segment was divided by the surface area of the segment, yielding a 
local wall flvix . This wall flux was then divided by the inlet flux. The procedure was 
repeated Avith 10 000 trials at the other tube end, Avith molecules initially at 295.0 K. 

Note that the wall flux for incoming 77. 4 K molecules beyond the tube center L/2 (where 
the temperature break occurs), going generally in the direction T^ — T 2 , is greater 
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Figure 11. - Variation of wall flux (normalized by flux entering tube) with axial position in 
tube. Tube length to radius ratio, 14.0; number of Monte Carlo trials, 10 000; gas to 
surface atomic mass number ratio, 0.21; step function wall temperature distribution with 
temperature change at midpoint. 


than the wall flux for incoming 295.0 K molecules beyond the tube center, going gener- 
ally in the opposite direction, Tg — T^. Thus, the difference in the average collision 
numbers KQ-yy(l, 2) and Kq^( 2,1) arises from gas-wall collisions after the molecules 
cross the wall -temperature break, which was placed near the center of the tube to ap- 
proximate experimental conditions. 

Figure 11 also shows that the wall flux is very nearly linear in the axial coordinate 
when Fq = 1.0 and all gas-wall reflections are diffuse. This well-known result was 
originally an assumption in Clausing’s integral equation approach to tube transmission 
probabilities and has been verified by many analytical studies since then (e.g. , refs. 15 
and 41). Also, the normalized flux at the midpoint of the tube when F^, = 1.0 is very 
nearly 0. 5, which also agrees with previous studies (e.g. , ref. 41) . 

Note in figure 11 that the two wall flux distributions when F^, is only 0. 3 are still 
fairly linear. Lobular reflections (forward scattering processes) tend to send incoming 
molecules down the tube more so than diffuse reflections, which results in smaller im- 
pact densities near the tube entrance plane (the left tube end in fig. 11) when F^, = 0.3. 
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Comments on Statistical Accuracy 


Throughout this study, because of the length of computer time required for calcula- 
tions involving lobular reflections in long tubes, a compromise had to be made between 
the number of cases studied and the accuracies of transmission probabilities for a given 
set of parameters. The 95-percent confidence intervals for are in the range 
e N = 0.007 to 0.0096 when K Q = 10 000 and are always small enough to statistically 
insure that > Qj 2* ^s “mentioned before, to halve the e-^, we must quadruple the 
number of Monte Carlo trials K Q or equivalently, quadruple the computer time. Dou- 
bling the number of trials to K q = 20 000 had little effect on answers in typical situa- 
tions, as seen by comparing cases 1 and 2 in table HI. 


TABLE III. - EFFECTS OF PARAMETRIC VARIATIONS IN PRESENT MONTE CARLO TUBE MODEL 


[Probability of cosine law reflection, 0. 3; step function wall temperature (see fig. 5(b)) unless otherwise noted.] 
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35 

.3340 

.3630 

1.087 

16.5 

15.2 

9 

.21 





6 

.0 





. 5335 

. 5390 

1.010 

7. 5 

6.8 

10 

.21 





20 

.0 

15 000 



.2617 

.2910 

1. 112 

24.0 

22.2 

11 

.21 



600.0 

14 

.0 

10 000 



.3130 

.3638 

1.162 

17.2 

15.7 


a All lobular reflections at 0 only. 
, max 

Partial ramp function (fig. 5(c)). 
c Fuli ramp function (fig. 5(d)). 


A trial case was investigated in which each lobularly reflected molecule was as- 
signed a reflection angle corresponding to the direction of the lobe pattern maximum 
6 v . For L/r = 14.0, = 0.3, a = 0.35, and p = 0.21, we found a T = 1.241 when 

all lobular emissions were at 0 max only. This number is different from a T = 1. 119, 
which resulted when sampling from the LKS distribution was done (see cases 1 and 3 in 
table III). Hence, although the computer time was less by about a third, the shorter 
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procedure of sending all lobularly reflected molecules off in the direction of the lobe 
maximum was not proper. Molecules reflected above and below 0 max apparently play 
an important role in determining the magnitude of a,p. Therefore, the more time con- 
suming process of sampling from the LKS distribution must be done. 


Dependence on Wall-Temperature Distribution 

The transmission probabilities and depend on the axial position of the 

step function wall -temperature break in the model as shown in table IV and figure 12 . 
Both Qjg and Q21 increase as the proportion of the tube that is cold increases, with 
Q21 being always larger than Q^- As the liquid level increases, a smaller part of the 
tube is hot. Hence, we would expect fewer cold gas - hot wall collisions. These kind 

TABLE IV. - DEPENDENCE OF TUBE TRANSMISSION 
PROBABILITIES ON POSITION OF WALL- 
TEMPERATURE DISCONTINUITY (STEP 
FUNCTION) CALCULATED BY USE 
OF PRESENT MODEL) 

[Tube length to radius ratio, 14.0; number of Monte 
Carlo trials, 10 000; temperature accommodation 
coefficient, 0.35; gas-to-surface atomic mass 
ratio, 0.21; probability of cosine law reflec- 
tion, 0. 3; Tj = 77. 4 K, T g = 295. 0 K.] 


Position of wall- 

Tube transmission 

Transmission 

temperature 

probability 

probability 

discontinuity, 

Q 12 

Q 21 

ratio, 

from tube end 

a T = Q 21 /Q 12 

a 0 

0.3008 

0. 3384 

1.125 

L 

.3138 

.3441 

1.097 

4 




L 

.3229 

.3614 

1.119 

2 




3L 

.3378 

.3677 

1.089 

4 




b L 

.3384 

.3745 

1.107 


a Entire tube at Tg = 295.0 K. 
^Entire tube at T 1 = 77. 4 K. 
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Liquid-nitrogen level height, fraction of tube length L 

Figure 12. - Variation of Q 12 and Q 21 with liquid-nitrogen level 
height using the present model. Tube length to radius ratio, 14.0; 
number of Monte Carlo trials, 10 000; temperature accomodation 
coefficient, 0.35; gas to surface atomic mass number ratio, 0.21; 
probability of cosine law reflection, 0.3. A step function, between 
Tj = 77.4 K and T 2 = 295.0 K, is used to simulate the nitrogen 
level. 

of collisions, considering lobular reflections, work to hinder cold gas motion toward the 
hot end. Therefore, Q.^ increases. We would also expect more hot gas - cold wall 
collisions as the level rises, for gas entering the hot tube end, which works to aid hot- 
gas motion toward the cold end. Hence, also increases. However, both values 

increase in such a way that a T = Q 21^12 rema i ns relatively unchanged. This is con- 
sistent with observed two-bulb behavior. 

A comparison of these values of arp (fig. 12) with some experimental values a^ is 
shown in figure 7. The theoretical values (for L/r = 14.0) are in the range of experi- 
mental values for the two-bulb data with L/r = 25.3. A detailed quantitative compari- 
son is not intended: calculations with L/r = 25.3 would require much more computer 
time. (However, additional results indicate changing L/r from 14. 0 to 20. 0 affects 
a^, only slightly . ) The difference between supposedly similar tubing samples 1 and 2 
in figure 7 points up the facts that knowledge of the nature of glass surfaces is incom- 
plete and that reproducibility of glass surfaces could be difficult to achieve experimen- 
tally. One could explain such a difference by assuming a larger F^, for sample 2, 
which might be expected if the inner surface of sample 2 was microscopically rougher 
than that of sample 1. 
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Both experimental results and the results of the model calculations indicate that 
even when the liquid -nitrogen level is at z = 0 or z = L and the entire tube is assumed 
to be at 295.0 or 77. 4 K, respectively, values of a E (or a,p) can still be greater than 
one. From the standpoint of thermal transpiration ratios, this implies that there can be 
a difference between two chambers joined by an aperture (which leads to a E = 1.00) and 
two chambers joined by a tube at the temperature of one of the chambers. In the model, 
even if the tube is at a single temperature, say, 295.0 K, that portion of the incoming 
77. 4 K gas that is reflected lobular ly can experience cold gas - hot wall reflections, 
which hinders passage through the tube. The 295.0 K gas entering from the other end 
can only have lobular reflections with T^ = 295.0/295.0 = 1.00. 

Another indication of the insensitivity of a^ to the tube wall -temperature distribu- 
tion between the 77.4 and 295.0 K chambers is the following: Instead of a step function 
for T g (z), as shown in figure 5(b), which led to arp = 1.119 (case 1, table in), we used 
the partial ramp function in figure 5(c). This led to a^ = 1.113 (case 4, table in). The 
use of a full ramp function, with T g (z) linear between z = 0 and L, led to a,p = 1. 127 
(case 5, table in). One might expect the step function to approximate the wall- 
temperature distribution with a heater wrapped around the tube just above the liquid- 
nitrogen level; and the partial ramp (or full ramp) function might better approximate the 
wall -temperature distribution without the heater. Although Edmonds and Hobson did not 
present heater - no heater results from two-bulb apparatus, they found in their elbow 
experiments that the heater either lowered the value of a E or left it unchanged. Our 
results indicate almost no change. 

Considering the calculated magnitudes of a^ and the insensitivity of a,p to tem- 
perature variations along the tube, the present model appears to simulate two -bulb data 
better than elbow data. 


Effects Involving the Nature of the Gas 

Hobson, Edmonds, and Verreault did most of their studies with elbow apparatus 
rather than two-bulb apparatus, since they were interested in obtaining and investigating 
the larger deviations from a E = 1.00. Bearing in mind the possibility of effects due to 
nonMaxwellian tube entry conditions, we will next consider several parameters which 
were experimentally investigated only in non -two-bulb apparatus. Unless otherwise 
stated, the subsequent calculations will be with F^, = 0.3, L/r = 14.0, p = 0.21, 

K q = 10 000, and the step function wall-temperature distribution in figure 5(b). The 
value F^, = 0.3 was chosen since it yielded a T values comparable with those obtained 
experimentally (fig. 7). 

Varying the species of the gas molecules can be considered in the model in three 
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I 


II 


ways: varying a q and/or the temperature accommodation coefficients; varying 

H = nig/m g ; and varying F^,, the probability of a diffuse reflection. As mentioned pre- 
viously, experimental results indicate that, as the mass of noble gases increases, mea- 
sured values of a E more closely approach 1.00. It has also been observed (e.g. , 
ref. 49) that thermal accommodation coefficients of gases on most surfaces generally 
increase as the gas mass increases. The effect of mass changes through accommodation 
coefficient changes was investigated by calculating Q 12 , Q 21 , and a T with = a L H 
a =1.00 (retaining p = 0. 21) and comparing with previous answers using ag = = 

a = 0.35. When a = 0. 35, we find a T = 1. 119, and, if a = 1.00, then a T = 1.137. 

(See cases 1 and 6, table III.) Increasing a to 1.00 has little effect on Q 12 , Q 21 , or 
a T . Other accommodation coefficient changes also result in relatively small changes in 
a T . For example, if a c = 1.00 and o? L = 0.35, then a T = 1.102 (case 7, table III). 
Thus, varying the two temperature accommodation coefficients for diffusely and lobular ly 
reflected molecules between 0.35 (an experimental value for helium on glass) and 1.00 
(the commonly assumed value for most gases on most surfaces) results in, at most, a 
2- to 3 -percent change in a T - 

To examine the effect of the mass ratio within the LKS distribution, values of Q^ 2 , 
Q 2 ^, and a T were calculated for p =0.33 and 0.21. Changing p from 0.21 to 0.33 
changed a T from 1. 119 to 1.087, a slight change in the correct experimental direction. 
(See cases 1 and 8, table in.) It was impossible to use a larger mass ratio appropriate 
for neon, for example. Using the assumption about the typical mass of a glass surface 
atom, the ratio would then be (20. 18/19. 2) = 1. 05. The LKS theory is valid only for 
p < 1/3, hence the limiting value p =0.33 was used. 

There is another mass effect which is probably more important than the direct effect 
of p in the LKS distribution in explaining (within the model) why a E decreases as m^ 
increases. It is plausible to expect that the probability of a diffuse reflection F^ in- 
creases with increasing m g . The more diffuse gas-surface reflections there are (i.e. , 
the larger F^, is) the less effect differential lobular scattering has on tube transmission 
probabilities. The scattering patterns observed by O’Keefe and Palmer (ref. 29) show 
that argon is reflected from polished glass surfaces much more diffusely than is helium. 
Generally, the larger the gas molecular mass, the more efficient is the energy transfer 
in a gas-surface interaction and the more likely a multiple collision or temporary ad- 
sorption is, which is usually thought to be followed by diffuse re-emission. However, 
current knowledge of gas-surface interactions is not complete enough to allow a reliable 
correlation of F^-, with gas mass. Varying F^, in such a way as to obtain a specific 
observed mass trend in the thermal transportation ratios would be of doubtful value with- 
out some theoretical foundation. 
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Effects Involving the Nature of the Surface 


Surface parameters also affect TTR values. A quantitative explanation for surface 
effects requires a more sophisticated tube model. Specifically, what is probably needed 
is a detailed correlation of F^, with gas and surface properties. When tubes are 
smooth, polished metal, Lund and Berman (ref. 48) found that transmission probabili- 
ties are fairly accurately predicted by theory and are not very temperature dependent. 
Consistent with this finding, Hobson (ref. 11) found for metal tubes that normalized 
ratios a^ in TTR experiments varied from unity by no more than several percent. At- 
tempting to simulate such situations within the LKS distribution, preliminary calcula- 
tions (meaning, with low values of K Q , the number of Monte Carlo trials) were done 
using a mass ratio p = 0.07, appropriate for helium on iron, nickel, or stainless steel. 
These preliminary calculations indicate as much or more difference in an< * ®21 

when p = 0.07 as when p = 0. 21. However, it is not certain whether the LKS model 
adequately predicts scattering patterns for p values much below about 0.10. Theo- 
retical lobes for p =0.02 (helium on platinum) are considerably thinner (narrower) than 
corresponding experimental lobes (fig. 4). 

It is likely that metal tubes, like increased gas molecular mass, lead to increased 
probabilities of diffuse reflection. This result might be expected by virtue of a deepened 
interaction potential between gas and surface molecules, or as a result of reflections 
within microscopic surface asperities. An increase in is also the obvious explana- 
tion for Hobson’s observations that in leached Pyrex tubing, Knudsen limiting law be- 
havior, Pj/Pg = (T^/Tg)* , was restored (ref. 10). Leaching refers to the dissolving 

action of HC1 on the inner surface of the Pyrex tubing, leaving a rough, porous surface 
layer. Not only should nearly complete diffuse reflection be insured by such surface 
roughness, but also virtually complete thermal accommodation, since each gas molecule 
can make many collisions within the porous layer before returning to the gas phase. 
Figure 10 shows that as F^, — 1. 0, "“Q21’ anc * a q< — 1-000- 

On the other hand, the presence of preadsorbed layers of Argon, Krypton, and Xenon 
on smooth Pyrex tubes had no effect on a^ for helium (ref. 11). In terms of the simple 
model, no detailed explanation can be offered now. Note that although the molecular 
weights of these preadsorbed gas layers are similar to those of many metals, the layers 
apparently behave like Pyrex surfaces instead of metal surfaces, in the thermal trans- 
piration effect. 


Effect of Tube Dimensions 

Another parameter found experimentally to affect TTR deviations is the tube length 
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to radius ratio. Edmonds and Hobson suggest that equation (1) is more closely ap- 
proached in tubes as the tube diameter is reduced (ref. 9) or as L/r is increased. For 
helium in a Pyrex tube of L/r =* 429, a value of a E = (Pi/P 2 )/R m = (0. 530/0. 512) ^ 
1.035 was found, substantially less than values for lower L/r tubes (ref. 9). How- 
ever, it was also found that a E = 1.00 when L/r is decreased to zero, or, when the 
tube becomes an aperture. Using the model presented herein with L/r = 6.0 yields 
a T = 1.010, indicating that such a tube is aperture-like (case 9, table III). If L/r = 

20.0, then a,p = 1.112, essentially unchanged from the value of a,p = 1.119 for L/r = 

14. 0, although both Q 12 and Q 21 were significantly reduced going to the higher L/r 
(case 10, table HI). Unf ortunately , because of excessive computer time, we were not 
able to investigate some of the capillaries of high L/r, say, 50 or 100 or more, used 
by Hobson and his co-workers. 

However, the effect of such long tubes in the model can be estimated after first dis- 
cussing short, aperture -like tubes. To obtain the TTR deviations in the model, it is 
important to have as nay hot gas - cold wall and cold gas - hot wall collisions as pos- 
sible; that is, trajectories of the type illustrated in figure 13(a), as opposed to those il- 
lustrated in figure 13(b) . Trajectories that do not lead to a difference in transmission 
probabilities (especially those in the first sketch of fig. 13(b)) could be expected to pre- 
dominate in tubes of small L/r, causing them to act as apertures in the thermal trans- 
piration effect. 

In very long tubes, not as many incoming gas molecules reach the vicinity of the 
wall -temperature change in the middle of the tube, even though the percentage of diffuse 



(a) Trajectories leading to difference in directional 
transmission probabilities. 
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(b) Trajectories not leading to a difference in directional transmission probabilities. 

Figure 13. - Sample molecular trajectories in tubes with step function wall temperature profiles shown. 





reflections is only 30 percent, as is assumed in most of the present calculations. In a 
very long tube, the number of diffuse reflections increases, which has the effect of 
erasing molecules' past history. So, even though hot -gas molecules are aided slightly 
near the middle of the tube and cold molecules are slightly hindered there, the L/r can 
be so large that the net separation effect on molecules successfully traversing the entire 
tube is lost, resulting in — Q21‘ ** ence > we believe the model would predict that 
a.rp — 1.00 as L/r — °°, assuming the wall -temperature discontinuity is near the middle 
of the tube . 


Effect of Incoming Hot Gas - Hot Tube Section Temperature 

The model indicates that the gas-surface temperature ratio is important in deter- 
mining the observed TTR deviations. To our knowledge, a thorough experimental study 
of the effect of varying gas and surface temperatures from the conveniently obtainable 
liquid nitrogen and room temperatures has not been made. In the model, Tg was 
changed from 295.0 to 600 K, resulting in a change in a,p from 1. 119 to 1.162 (cases 1 
and 2, table III). This trend is the same as indicated in Hobsons' three-stage accommo- 
dation pump results using Tg = 600 K instead of 295.0 K (ref. 11). As previously 
stated, our single -tube results cannot reliable be extrapolated to predict pump ratios 
across a three -stage pump because of uncontrolled and unpredictable variations in sur- 
face conditions (ref. 11), and possible nonMaxwellian tube entry conditions in such a 
pump. 


Predicted Axial Variation of Gas Temperature 

The technique for calculating transmission probabilities readily allows calculation 
of the average temperature of incident gas molecules at any axial position, as well as 
the average temperatures of the molecules leaving the tube through the two end planes. 

It is important to note again that we are ascribing to each gas molecule the property of 
temperature, regarding temperature as a measure of average gas molecular kinetic 
energy. Hence, the following results should be regarded only as first-order indications 
of how gas temperature would vary axially, with partial LKS lobular scattering and low 
temper ature accommodation at the tube wall. 

Shown in figure 14 are the average temperatures of gas molecules striking the tube 
wall as functions of axial position. The temperatures are calculated for 10 000 mole- 
cules entering the cold end (fig. 14(a)) and then for 10 000 molecules entering the hot end 
(fig. 14(b)), with F^, = 0.3 and a = 0.35. Then, these calculations were repeated with 
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on 



(b) All gas molecules started initially at the T 2 = 295.0 K tube end . 


Figure 14. - Average temperature of gas molecules striking the tube wall as function of 
axial position. Tube length to radius ratio, 14.0; number of Monte Carlo trials, 

10 000 (in each direction); gas to surface atomic mass number ratio, 0.21; wall tem- 
perature step function at L/2 between and T 2 tube sections. 


F^, = 1.0 and a = l.O. When = 0.3 and a = O. 35, the originally hot gas has only 
cooled down to about 110 K when it strikes the walls near the cold end (fig. 14(b)), and 
the originally cold gas has only heated up to about 270 K near the walls of the hot end 
(fig. 14(a)). 

The average temperatures of molecules leaving both tube ends, summarized in 
table V, reveal even more discrepancy between gas and wall temperature: When 


34 



TABLE V. - COMPARISON OF CALCULATED AVERAGE TEMPERATURES OF GAS 


MOLECULES LEAVING TUBE THROUGH EXIT AND ENTRANCE PLANES 


[Tube length to radius ratio, 14.0; assumed helium to Pyrex atomic mass ratio, 
0. 21; number of Monte Carlo trials at each end, 10 000. ] 


Condi- 

tion 

Assumed 
probability 
of cosine 
law 

reflection, 

F c 

Assumed 

temperature 

accommodation 

coefficient, 

a 

Molecules 
entering at 
tube end 

(a) 

Calculated 
average 
temperature 
of molecules 
leaving through 
exit plane, 

t ex> 

K 

Calculated 
average 
temperature 
of molecules 
leaving through 
entrance plane, 

t en> 

K 

1 

1.0 

1.00 

T 1 

273.8 

78.8 

2 

1.0 

1.00 

T 2 

96.9 

293.7 

3 

0.3 

0.35 

T 1 

248.2 

87.6 

4 

.3 

.35 

T 2 

141.3 

289.1 


a T x = 77. 4 K; T 2 = 295.0 K. 


Fq =0.3 and a = 0. 35, the average temperature of the orginally hot gas which even- 
tually exits the cold end is about 143. 8 K, while the average temperature of the originally 
cold gas which eventually exits the hot end is about 246. 4 K. Even when F c = 1. 0 and 
a =1.0, these exit temperatures are 96.9 and 273. 8 K, respectively. 

This difference between the gas and the wall temperature could be important when 
trying to assess the behavior of the gas in the chambers (bulb or elbow) at the tube ends. 


Effects Peculiar to Non-Two-Bulb Data 

The model we have discussed assumes that the gases in the reservoirs or the tube 
ends are Maxwellian. This assumption yields the cosine distribution at the entry planes 

and the ^2-nm^.T^ entry rate, which leads to equation (2). When the reservoirs 

are large bulbs, one would expect the Maxwellian assumption to be valid. However, if 
the reservoir is an elbow-shaped, ill-defined tubular region, the gas entering the tube 
of concern could take on a nonMaxwe Ilian character. This may be the origin of the be- 
havior peculiar to elbow apparatus in Hobson's experiments: larger values of a E (in 
the approximate range 1.08 to 1.31 for helium in Pyrex tubing) and a sensitivity of a E 
to liquid -nitrogen height up the tube. 
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If an elbow-like reservior is assumed, it is not at all certain the the directions of 
the molecules on the entrance plane are described by a cosine distribution. Because of 
a lack of detailed information about the geometric shape and surface properties of the 
elbow joint beyond the entrance plane, it is difficult to estimate the directions of enter- 
ing molecules. Dayton points out that in calculating tube conductances, one must con- 
sider the beaming pattern across the entrance as well as the temperature and back- 
scattering characteristics of surfaces beyond the exit of the tube (ref. 50). It is well 
known that unlike the cosine efflux pattern from a small hole in a bulb reservoir, tubes 
even with an L/r as small as 2. 0 exert a considerable beaming effect in an axial direc- 
tion on exiting molecules, as shown, for example, by Chubb (ref. 51). The only work on 
rarefied flow through elbows (by Davis, ref. 16) deals only with transmission probabili- 
ties, not efflux patterns, but one could reasonably expect some type of beaming effect 
here too. 

In the present model, the use of two different entry distributions designed to beam 
incoming molecules down the tube more axially than a cosine distribution on the entrance 
plane did not result in the larger values of a T suggested by elbow experiments. 

Another theoretical difficulty inherent in the elbow reservoir is that it may not allow 
enough gas-surface interactions to insure that the gas temperature equals the wall tem- 
perature. Results in the previous section indicate that the gas temperature can be sub- 
stantially different from the wall temperature near the tube ends, for gases with small 
accommodation coefficients and for a tube L/r of 14.0. Consequently, incorrectly as- 
sumed temperatures in equation (2), as well as a difference in and Qgp help de- 

termine the behavior of the deviations of P^/Pg from the Knudsen limiting law (eq. (1)) 
observed in the elbow data . 

At present, however, not enough is known about nonMaxwellian entry conditions to 
warrant detailed modifications of the current model, in an attempt to explain the elbow 
TTR effects. 


Comments on Accommodation Pumping 

The action of the accommodation pump of Hobson (fig. 6(c)) can be qualitatively ex- 
plained by considering transmission probabilities in the present model. Note that the 
pump is constructed of alternating smooth and rough (leached) sections of Pyrex tubing. 
Our calculations show that the transmission probabilities of molecules moving in the hot 
(top) to cold (bottom) direction in the smooth tube sections are greater than those found 
for molecules moving in the cold-to-hot direction because of temperature -dependent 
lobular scattering. On the other hand, the rough, leached tube sections, probably having 
an ^ 1.0, would not be expected to sustain much lobular scattering, and thus, in 
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these sections, transmission probabilities in the hot-to-cold direction would nearly equal 
those in the cold-to-hot direction. Thus, in effect, each adjacent pair of smooth and 
rough tubes will, because of the differential transmission probabilities, act as a pump. 

In the accommodation pump, each such stage acts in series with other stages to multiply 
the overall pumping effect. 


SUMMARY OF RESULTS 

A model of free molecular gas flow in tubes was developed that combines a lobular 
reflection distribution with conventional cosine law (diffuse) reflection at the wall using 
a Monte Carlo approach to the solution. At each gas-wall collision, an incoming mole- 
cule is either reflected diffusively with probability F^, or it is reflected in a direction 
chosen from the Logan-Keck-Stickney lobular reflection distribution, with probability 
1 - F^,. Assumed temperature accommodation coefficients determine the temperature, 
or average kinetic energy, of the gas molecule after the interaction, in either case. 

This lobular distribution is a function of incident gas and surface temperatures, inci- 
dence angle, and gas-to-surface atomic mass ratio and has successfully predicted a 
number of observed gas molecular -beam surface -scattering patterns. 

The model was used to explain previously observed deviations in the free molecular 
thermal transpiration ratio from the Knudsen limiting law. These deviations suggest 
that when two reservoirs at temperatures T^ and Tg are joined by glass tubing , gas 
molecules are more likely to traverse the tube in the hot-to-cold direction than in the 
cold-to-hot direction; that is, Qgj > when Tg > T^, where is the tube trans- 
mission probability from reservoir i to reservoir j. 

The major results from the calculation of and using the model are the 

following: 

1. Theoretically calculated values of a^, = Q 21 /Q 12 are the ran S e 1-09 to 1. 14 
for helium when F^ is assumed to be 0.3 (the same for all gas-wall collisions), when 
T^ = 77. 4 K and Tg = 295.0 K, and when the tube L/r = 14.0. This is consistent with 
some previous thermal transpiration ratio data obtained with Pyrex tubing joining two 
glass bulb reservoirs (two-bulb data) . 

2. When part of the length of the tube is assumed to be at T^ and the remainder at 
Tg, the average number of gas-wall collisions per incoming gas molecule is greater in 
the T^ — Tg direction than in the Tg — T^ direction. The difference arises from colli- 
sions after the molecules across the wall -temperature break, which was placed near the 
tube's center to approximate experimental conditions. This means that initially hot 
molecules make slightly fewer wall collisions, on the average, than initially cold mole- 
cules starting from the opposite tube end. 
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3. Values of arp are relatively insensitive to changes in the wall -temperature dis- 
tribution between the and Tg tube ends. Also, a^>1.00 even when the entire 
tube is at or Tg. Both findings agree with previous two -bulb data. 

4. Values of a^ are fairly insensitive to varying the thermal accommodation coef- 
ficients between 0.35 and 1.00. 

5. The model confirms the following experimental trends: The transmission proba- 
bility ratio arp 

(a) decreases to 1.00 as the tube length -to-radius ratio approaches that of an 

aperture, zero 

(b) decreases toward 1.00 as the tube wall becomes microscopically rougher 

(assuming F^, increases as roughness increases) 

(c) increases as the hot gas - hot tube section temperature increases 

(d) decreases as gas molecular mass increases (by changing the gas-to-surface 

mass ratio; but, more reasonably, by assuming F^, increases as m^ in- 
creases) . 

6. The model is qualitatively consistent with Hobson's accommodation pump studies, 
which used such tubes with transmission probability directional differences. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 19, 1971, 

113-31. 
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APPENDIX - SYMBOLS 


A 1> A 2 

a E 


B r B 2 

cos(A , cos(Bj) , cos(C.j) 


*L 

f(u) 

K GW^’^ 


K o 

k 

L 

LKS 

m g 

m s 

P A> P B 

P i 

p N (0 r ) 

P(0 r ) 


Q 

Q 

R 

R. 


D 


i] 


tube entrance planes 

experimental thermal transpiration ratio, (temperature 
normalized), (P 1 /P 2 )/R m 

transmission probability ratio, Q21/Q12 

quantities in LKS equation (eq. (3)) 

direction cosines of particle trajectories, with x, y, z axes, 
respectively 

probability of cosine law (diffuse) reflection 
probability of lobular reflection, 1 - 
probability density function 

average number of gas-wall collisions per incoming gas mole- 
cule at T. , moving generally in the T. direction; i,j = 1,2 

number of Monte Carlo trials 

Boltzmann constant, 1.38054x10“^ JK“* 

tube length, m 

L ogan -Ke ck -Stickney 

gas atomic mass number 

surface atomic mass number 

pressures at end of an accommodation pump 

gas pressure in region i; i = 1,2 

wp/Pi'W 

probability of reflection into d$ r at 6 V , given by LKS 
equation (eq. (3)) 

probability of falsity of inequality in eq. (10b) 
direct tube transmission probability 

tube transmission probability in i — j direction; i, j = 1,2 
thermal transpiration ratio (P^/Pg) 

random variable uniformly distributed on the interval (0,1) 
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(t 1 A 2 ) 1/2 

tube radius, m 

set of random numbers 

average temperature of gas molecules leaving tube through 
entrance plane, K 

average temperature of gas molecules leaving tube through 
exit plane, K 

gas temperature, K 

gas temperature in region i, K; i = 1,2 

temperature of gas molecules leaving and approaching tube 
wall, respectively, K 

surface -to-gas temperature ratio, T g /Tg 

surface temperature, K 

wall -temperature distribution as function of axial coor- 
dinate z, K 

thermal transpiration ratio 
variable in probability density function 
limiting values of u 

coordinates of particle position in tube model 
temperature accommodation coefficient 

temperature accommodation coefficient for cosine reflections 
temperature accommodation coefficient for lobular reflections 
Chevyshev confidence limit 

95-percent confidence limit for Q^. assuming normal distri- 
bution of Monte Carlo answers 

random variable 

angle of incidence, measured from surface normal, deg 
value of $ r at which P(# r ) is maximum 
angle of reflection, measured from surface normal, deg 
specular reflection direction, = 9 ., deg 
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